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A stream of telephone calls is submitted to a group of trunks, the first- 
choice group, according to a recurrent process. We allow balking on this 
trunk group; i.e., if a call finds k of the first-choice trunks busy it may be 
served, with probability pk , or may fail to be served, with probability gu ■ 
A call which fails to receive immediate service on the first-choice trunk group 
is submitted to a second-choice trunk group, the overflow group. We also 
allow balking on the overflow group. Calls which fail to receive immediate 
service on the overflow group are lost to the system. Holding times have 
negative-exponential distribution . 

We give methods for finding the joint distributions of numbers of busy 
trunks on the first-choice and overflow groups, at overflow instants (i.e., 
instants at which calls are submitted to the overflow group), at arrival in- 
stants, and at arbitrary instants. We consider the transient as well as the 
limiting distributions (and demonstrate the existence of the limiting distri- 
butions) . 

The methods developed are illustrated by several examples. Numerical 
results arc given for the blocking in (he particular case that the first-choice 
group constitutes a random slip, while the overflow group is full-access 
(common) . 

I. INTRODUCTION 

i.i Balking and Overflow Traffic 

A telephone call is submitted to a group of m trunks. This call may 
fail to occupy a trunk, even though not all m trunks are busy. There 
may be a number of reasons for such a failure, e.g. : the calling line may 
not have access to any idle trunks, some equipment other than the 

* This paper represents part of a doctoral dissertation submitted to the Sub- 
committee on Applied Mathematics, Columbia University. 
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trunk itself may be required to complete a connection and this equip- 
ment may be busy, or the m trunks may be merely first-stage links in 
a connecting network and there may be no free path through this net- 
work. Whatever the cause of the failure, we shall say that the submitted 
call balks (although the word is perhaps more appropriate in queueing 
theory applications). In this paper we shall restrict ourselves to the case 
in which the probability of balking depends only on the number of 
busy trunks: if an arriving call finds k trunks busy, it is served, with 
probability p k , or balks with probability q k (p k + q k = 1). If all trunks 
are busy, an arriving call cannot be served, and therefore q m = 1. Thus 
we subsume blocking under the term balking. 

The traffic which overflows from a trunk group with balking has 
different characteristics from that which overflows from a full-access 
group. [By a full-access trunk group we mean one for which q k = 
(k < m), q m = 1.] Suppose recurrent traffic is submitted to a full-access 
group (when we refer to recurrent input traffic we mean that the inter- 
vals between arriving calls are independent, identically distributed 
random variables). Suppose further that the holding times of calls have 
negative-exponential distribution. Then, as Conny Palm 1 has shown, 
the overflow traffic is also recurrent. This is not the case for traffic 
overflowing from a trunk group with balking. 

The traffic which balks on the first-choice group may be submitted 
to an overflow group of, say, M trunks. There may also be balking on 
the overflow group. Now L. Takacs 2 has treated in detail the process 
of numbers of busy trunks in a trunk group with balking to which a 
recurrent stream of calls of negative-exponential holding times is sub- 
mitted. Thus, if the first-choice group is full-access, we know how to 
describe what goes on on the overflow group. However, if the first-choice 
group is not full-access, the stream of calls submitted to the overflow 
group is not recurrent, and therefore further analysis is required to 
describe the process of numbers of busy trunks on the overflow group. 
We attempt to treat this problem in the present paper; in so doing, we 
are led to consider the joint distribution of numbers of busy trunks on 
the first-choice and overflow groups, which is also of interest in itself. 

1.2 Mathematical Description of the Problem, and Some Notation 

Calls are submitted to a group of m trunks, the first-choice group, at 
successive instants n , t 2 , • • • , t„ , • • • . The interarrival times, 6 n = 
Tn _ Tnl ( n = 2, 3, 4, • • • ), are independent, identically distributed 
random variables with common distribution function 

P{0 n ^ x] = F(x), 
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and we specify further that P\n <; x\ = F(x). We assume that the 
\6„\ are not lattice variables (i.e., that the interarrival times are not 
confined to multiples of a constant), that F(0) = and that 

< a < co, 

where 



a = f xdF(x) 



is the mean interarrival time. 

Note that the class of recurrent inputs just described includes, among 
others: Poisson arrivals, equally spaced arrivals, and, as previously 
remarked, arrivals which are themselves overflows from a full-access 
trunk group to which a Poisson process of calls with negative-exponential 
holding time is submitted. 

If the nth call receives service, then its holding time is a random varia- 
ble, Xn • The \xn] are independent and identically distributed, with 
common distribution function 

( 1 - e~ T for x > 
PiXn ^ x) = 

[0 for a < 

and are independent of the arrival process {t„}. 

Note that we are measuring time in units of the mean holding time; 
thus a = \/a is the submitted traffic in erlangs. 

An arriving call which finds k trunks of the first-choice group busy is 
served with probability p k , or balks with probability q k . We have 

Pk + Qk =1 (k = 0, 1, ••• , m) 

q m = 1. 

A call which balks on the first-choice group is immediately submitted 
to a second group of M trunks, the overflow group (we allow the case 
M = co ). We denote the sequence of instants at which calls are sub- 
mitted to the overflow group by { T N ) (N = 1, 2, 3, ■ • • ). If such a 
call finds K trunks of the overflow group busy, it is served, with prob- 
ability K , or balks, with probability II K . We have 

G K + H K : = 1 (K = 0, 1, ••• ,M) 

II M = 1 (if ^1/ < co). 

We make the following plausible restriction on the balking proba- 
bilities 
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p k > for k < m 
G K > for K < M. 

A call which balks on the overflow group is said to be blocked. It im- 
mediately disappears from the system and is not resubmitted; i.e., lost 
calls are cleared. 

We now define the following random variables: 

£(0 = number of busy trunks on first-choice group at time t 

& = £(r„ -) 

£ n ° = £(Tn — ) (the superscript "o" means "overflow".) 
S(0 = number of busy trunks on overflow group at time, t 

En = S(r„ -) 

t-t o tV rn \ 

AN — A{1 N )• 

We also define the following probabilities, which it will be our object 
to determine: 

P{6r - = h, 2/ = K) = P°(k,K,N) 

lim P°(k,K,N) = P°(k,K) 

N-*oo 

P{& = k, Sn = K\ = P(k,K,n) 
\imP(k,K,n) = P(k,K) 

n-»°o 

P{$(0 = k, SCO = K\ = r(KK,t) 

YunP(k,K,t) = P*(k,K). 



When one of the variables k or K in one of these probabilities is not 
written, it is understood to be summed over, e.g. 



P(k,t) = T,P(k,K,t). 
A quantity of particular interest in applications is the blocking 

We shall also have occasion to refer to the blocking on the first-choice 
group 

b = £ q k P(k). 
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Further notation will be introduced as it is needed. The notation will 
as far as possible conform to that of Takacs. 2 We shall, when possible, 
use lower-case letters to refer to the first-choice group and the corre- 
sponding capital letters for the overflow group. Equations of Ref. 2 
will be denoted by a T: e.g., "(T44)." We note here only the following 
definitions: 






C -Ur^L (ft-1) 

j-i 1 - <PU) 

cw-n . *'J; J 2n (c-w-i). 

1.3 Previous Results 

Let us denote the interoverflow times by ®/t = T N — T N ~\ . As we 
have mentioned, if the first-choice group is full-access, the j®^} are 
independent and identically distributed. In this case let us denote their 
common distribution function by 

G(x) = P{S N ^ x) 

with Laplace-Stieltjes transform 

7 (s) = I e-" r dG(x). 
Jo 

Takacs solves a recurrence of Palm to obtain 

±( m ) l 

y^ (m + 1 \ ! 



r / Cr-i(s) 

A. Descloux gives convenient recurrence formulas for calculating 
y(s) and the moments of G(x) in the case of Poisson input, i.e., when 

fl - e~ ax (x ^ 0) 
F(x) = 

[0 (x < 0) 

Some results exist for P(k,K) in the case of Poisson input [for which, 
and only for which, as we shall see, P*(k,K) = P(k,K)]. The first of 
these is due to L. Kosten. He considers a full-access first-choice group 
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and an infinite full-access overflow group. Let us denote binomial 
moments with respect to the overflow group by 

U(k,R) = Xj, (*£\ P(k,K). 
Then Kosten finds 

^^' ^g) - (2) 

(See also the appendix by J. Riordan to a paper of R. I. Wilkinson. 6 ) 
The polynomials in (2) are denned by 

k-i 

(3) 




so that CV'(a) = a/k\, if we agree that ( J - 1. J. Riordan (Ref. 

7, p. 120) remarks that these polynomials are closely related to the 
Poisson-Charlier polynomials C»(x,a); in fact 

C K \a) = C k {-R,a). 

E. Brockmeyer, 8 N. Bech, 9 and K. Lundkvist 10 consider a problem 
which differs from Kosten 's only in that M is finite (G M = 0). Brock- 
meyer finds 

P(k,K) = "f (-lVY B+K (* + *) C K+a "- s (a) (4) 



where 



Y s = E (-D y_S (i _ J) *' (S - 1, 2, • • ■ , M) 



Y»= * 



d m + M (a) 

aj = c^y ' cA^) £ v - 1) Co '" +L(a) - 

We do not consider here more complicated trunking situations (graded 
multiples, alternate routing arrangements in which the overflow group 
is at the same time the first-choice group for other sources of traffic). 
See, however, Wilkinson, 6 and R. Syski (Ref. 11, chapters 7, 8, 10). 

Takacs 2 gives, for arbitrary q k , methods of finding P(k,n), P(k), 
P(k,t), and P*(k). Thus in what follows we shall take the attitude that 
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everything we need concerning the first-choice group only is, in principle, 
known. 

i A An Example 

This paper grew out of the following problem, in which both balking 
and overflow are involved. Subscriber lines are connected to the in 
trunks of the first-choice group in such a way that each line has access 
to only 7 of them. We refer to a particular set of 7 trunks as the access 
pattern for a particular line or group of lines. Equal traffic is submitted 



to each of the I I possible access patterns. When a connection is made, 

any idle trunk in the subscriber's access pattern is equally likely to be 
selected. This arrangement is referred to as a random slip, or Erlang's 
ideal grade. It is easy to see that the balking probabilities are 



q k = 0, for ^ k < 7, and 
Qk = 7^4 , for 7 ^ k ^ m. 



Traffic which balks on the first-choice group is submitted to a full- 
access overflow group of ]\I trunks. If a call is blocked on the overflow 
group, it is lost. 

Such an arrangement may be economically desirable. The average 
traffic carried per trunk (for a given blocking probability, B) is less than 
for a full-access group of m + M trunks, but the traffic per crosspoint 
is greater. Knowing the costs of trunks and of crosspoints, and given 
m + M and the desired value of B, one wishes to select 7 and m so as 
to minimize the cost per unit of carried traffic. We shall give some nu- 
merical results for this arrangment. 

II. THE STATE OF THE SYSTEM AT OVERFLOW INSTANTS 

2.1 Transient Behaviour 

Unless the first-choice group is full-access, the overflow process { TV} 
is not recurrent and the sequence j Eat ) is not a Markov chain. However, 
the sequence of pairs of random variables {£*", Ea-°} is a homogeneous 
Markov chain. This may be seen as follows. Suppose we know that 
£(T N — ) = k and H( T N — ) = K. T N is an arrival instant; because the 
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arrival process is recurrent and independent of the holding times, the 
history of the system before T N has no effect on the epochs of future 
arrivals. T N is an overflow instant; whether or not the overflowing call 
is accepted by the overflow group depends only on the value of K. Be- 
cause of the exponential distribution of holding times, the stochastic 
behaviour of the system after T N is independent of the ages of calls in 
progress at T N . Thus the values of £(7V— ) and S( T N — ) determine the 
whole future stochastic behaviour of the system. Therefore we are led 
first to a consideration of the probabilities P°(k,K,N). 

If £(J) = k, S(0 = K, then we say that at time t the system is in 
the state (k,K). The values of £/ are limited to those k for which q k > 0. 
We denote the set of such integers k by a. As initial conditions we take 
£(0+) = i, 2(0+) = I < co. (It is not required that i £ a.) Under 
these initial conditions, we seek P°(k,K,N) for fc 6 G; K - 0, 1, 2, • ■ ■ ; 

N =1,2,3, ••-. 

Let us now define the following quantities : 

G ik (x) = PlH„ +l ° = k, 0„ +1 ^ x | £(T„+) = j\ 
= -PiW = K ©*+i ^ s|Sjv = j\ 
= PUi° = fc, T x £z|*(0+) =iJ 
with Laplace-Stieltjes transform 

y jk (s) = f e" x dG jk (x) 
Jo 

U°(k,R,N) = Z h (£\ r°(k,K,N) (R = 0, 1, • • • , M) 

F°(A-,fl,A r ) = ]£ (f) G K P°(k,K,N) (R = 0, 1, ■ • • , M) 

F°(fc,-1,A) - 0. 

We may now state: 

Theorem 1: The distribution P"(k,K,N) is uniquely determined by 
the binomial moments U°(k,R,N); the latter are determined by 

u°(k,R,i) = (jj.)*y*(fl) ( 5 ) 

U°(k,R,N + 1) - ETi*(fi)TO^) + ^°0',« - W]- (6) 
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Proof: The transition probabilities for the homogeneous Markov 
chain { £ /, E/j are given by 

p°(j,J; k,K) = P{ W = k, H„ +1 ° = K | & = j, -E N ° = J\ 

= I P (S.v-h = K | Zh = J, A '+i = x) dG jk (x). 
jo 

It is easy to see that 

P\ Aff+i" = K | Sat" = J, ©at +1 = x\ 

J "i 1\ -J-K/1 -x\J+l-/r 



= #/ ( t^ ) c ' ( 1 - e 



Thus 



hff/(i) «"**(i - e - j 



p°(j, J; k, K) = jT <flfe(z) fa/" 7 J l ) e-* K (l - e~* 



~ HAi)e- xK (] - e~ z ) 



—I \ J+l-K 

(7) 



Now 

m .1/ 

P°(k,K,N + 1) = E Y,V(hJJ<,K)P°(j,J,N). (8) 



y=o y-o 



Substituting (7) in (8), and taking the Rth binomial moment with 
respect to the overflow group, we obtain 

U°(k, R, N + 1) - £ £ f dGj»(aO ftf/" 7 + lN ) 
+ ff/(^)] *""** P°(j, J, N) 

- £. S ** ( * > [© + Gj (r - i)] P °°' J > N) 

= E t*(«) Itf'O', K, #) + V°0",ff - 1, N)], 



Jta 



which is (()). 
For N = 1, wc have 
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P°(k,K, 1) = [ dG ik {x) (£) e- xK (l - e-*) 1 -* 



so tliat 



(9) 



U°(k, R, 1) = jT dG ik (x) Q •* - Q Y*(i2), 

which is (5). 

From the definition of U°(k,R,N), we have 

£ (-i) K (£)Er(* f JB,tf) 

-£<- i >**S)iffi r * / ' ir> ' 

Now, for any finite N the double series on the right contains a finite 
number of terms, even if M = °o . This is so because 

P"(k,J,N) = for k + J ^ i + / + N, 

and we have assumed / < oo . 

Thus the double series can be rearranged, and one obtains readily 
that the binomial moments determine the probabilities according to 

P°(k,K,N) = jt (-I)*'* (*£)u°{k,R,N). (10) 

In (5) and (6), the quantities y,k(R) occur as coefficients. We regard 
these coefficients as known because they can be expressed in terms of 
certain quantities determined by Takacs. Let 

M ik {x) = E {number of t„ in (o,x] for which £„ = k \ £(0+) = i], 

with Laplace-Stieltjes transform 

Jo 

Takacs gives a method for finding the m«(s) [(T70), in which, however, 
the index i is implicit]. The way in which the quantities n Jk (R) deter- 
mine the y jk (R) is expressed in the following lemma (in which, it is to 
be noted, values of the indices j,k, etc. are no longer restricted to the 

set a). 

Lemma 1: Define M ik °(x) = E {number of T N in (0,x] for which £/ = 



OVERFLOW TRAFFIC 55 

k | 5(0+) = i), with Laplace- Stieltjes transform 

«/(•) = f e~ 8x dM ik "(x). 

Let n°' be the square matrix with elements /i/*°(/2); j,k = 0, 1, • • • , m. 
Let y R be the square matrix with elements y ik (R);j,k = 0, 1, • • • , m. 
Then, for R = 1,2, ••• , 

7 = n (L + m ' ) (11) 

where E is the (w + 1) by (m + 1 ) unit matrix. 
Since, obviously 

n Jk °(R) = q k n Jk (R), (12) 

(11) provides the desired relation between the yj k (R) and the fij k (R). 
Proof: We shall first show that 

HI 

*»(R) = yAR) + ZtM^WW (13) 

/=0 

for R = 1, 2, ••■ . 

Suppose £(0+) = j, and consider a given /?-tuple of trunks on the 
overflow group which are all busy at / = 0+. If 7\ = x, the probability 
that the overflow at 7\ will find this 22-tuple still busy is e~ Rx . 

Thus 



y Jk (R) = [ e~ Rx dG jk (x) 

Jo 



is the probability that this 7?-tuple is still busy at T x and that £(!Ti — ) 
= k. 

Again, if this 72-tuple remains busy just until t = x, the expected 
number of overflows from k to find it busy is M jk °(x). Therefore the 
unconditional expectation of the number of overflows from k to find 
it busy is 

f M jk °(x) d(l - e~ Rl ) = f e- RT dMjt'(x) = W *°(ff). 

•'o •'0 

Denote (temporarily) by \tij k °(R) | /] the expected number of overflows 
from k to find this i^-tuple still busy, on the condition that £{1\ — ) = I 
and the 72-tuple is still busy at t = T x — . 

Then, by the principle of total expectation, 

/**•(«) = E fo» # («) I ihdRh (14) 
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Now because of the exponential holding-time distribution 

M(B) I *] = Mu-°(tf) for I 9* k (15) 

and 

[n Jk °(R)\k] = l +m°(R). (16) 

Substituting (15) and (16) into (14), we obtain (13). Equation (13) 
may be written 

H ' = 7 + 7 M • ( 17 ^ 

Thus, to prove the lemma, it remains to show that {E + /» ' ) IS 
nonsingular. 
From (17) 

(E - y )n = 7 • 
Therefore 

(tf-7 - ) •(* + *'•*) = # 
det (£-7*) -det(JS? + M° ,fl ) = I- 

Since clearly both det (E - 7") and det (E + //'*) arc finite (for 
li > 0), it follows that det (# - 7") * and det (E + m°'*) ^ <>, 
which completes the proof of the lemma. 

We note, for later use, that we have also shown that 

M°' K = (E-y R )"Y. (18) 

We need a separate method for finding 7/*(0), the above argument 
being invalid because Mja-°(0) = <=° for all k £ ft. 

We notice that 7>a-(0) = (?,*(«) = P{£(7\-) - fc|S(0+) = /}. 

The quantities 7y*(0) are determined by the following system of 
equations: 



(19) 



7i*(0) = q* [ dF(x) 0j\ <f ta (l - e^Y'" + g p, tWO)- 

• jT dF(x) A') <r ,x (l - e-)^' 0', ft - 0, 1, • • ■ , to). 

This may be seen as follows: 

The event {£( T%— ) = k] can occur in these mutually exclusive ways: 
(i) the first arrival after t = encounters k busy trunks on the first- 
choice group, with probability 
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[ dF(x) K)«- te (l - e 

Jo \/b'/ 



\j-k 



and overflows, with probability qk ; 

(ii) the first arrival after t = encounters I busy trunks and does 
not overflow [so that £(7\+) = Z + 1]; the next overflow following 
this occurrence is from k [probability 7i+i,fc(0)]. 

For each k, (19) is a set of linear equations in the 7#(0). These 
equations determine the 7i*(0) uniquely if the coefficient matrix is 
nonsingular (for each k). Call this matrix A (k) . If we can show that 
I An*' \ > ^2 Aji' k) for each j, it will follow from the theorem of Levy- 

Hadamard-( Jerschgorin (lief. 12, p. 79) that det A m ^ 0. That is, 
we want to show that 

£ Pl f M dF(x) ({) e- lx (l - e-y- 1 < 1. (20) 

The left side of (20) is evidently strictly less than 

£ ( <lF(x) KJ e- lx (l - e~ x y- 1 = 1, for each j, Q.E.D. 

Equations (5) and (0) may he solved, in some cases, by means of 
generating functions. 
Let 

U°(k,R,w) = £ U°(k,R,N)w N 

.V=l 

V°(k,R,w) = £ V°(k,R,N)w N 

A = l 

Note that it follows from (10) that 

£ P°(k,K 1 N)w N - £ (-1)*"* (^)u°(k,R,w). (21) 

AT=1 H=K \A/ 

From (5) and (6) we obtain 
U°(k,R,w) = wi(^y ik (R) + ZyAR)[U°(j,R,w) 

y (22) 

+ V°(j,R - 1,«)] 
We illustrate the use of (22) by a simple example. 
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Example 1: 

If the first-choice group is full-access (the only element of ft is m), 
then U°(k,R,N) and V (k,R,N) vanish except for ft = m. For sim- 
plicity, we assume that i = m; then the only relevant element of the 
matrix 7" is y,nm(R), and (22) becomes: 

U°(m,R,w) m wy mm (R)UV\ + U°(m,R,w) + V°(m,R- l,w)j, 
whence 

^ * •> = T^Ek [© + y " (M ' * " '• w) ] ' (23) 

7mm(s) is * ne Laplace-Rtieltjes transform of the interoverflow-time 
distribution, i.e., it is just the function 7 (s) given by (1). Thus (23) 
is exactly equivalent to (T32), and merely serves to illustrate our remark 
(Section 1.1) that if the first-choice group is full-access, we can use the 
methods of lief. 2 to describe the behaviour of the sequence \En°). 

2.2 The Limiting Distribution P°(k,K) 

Theorem 2: The quantities P°(k,K) = Urn P°(k,K,N) exist, are 

JV-.W 

strictly positive, form a probability distribution independent of the initial 
slate, and are uniquely determined by the binomial moments U°(k,R) = 

E( p )P°(k,K); the latter are determined by 
k=r\R/ 

U°(k,R) = q k E nAR)V°(j,R - 1) (R - 1, 2, • • • , M) (24) 



and 



where 



U°(k,0)= q -^- (25) 



V°(k,R) = Jg(|) G K P°(k,K). 



Proof: We first show the existence of the limiting distribution. 

In this section, we use theorems given in Feller, 13 chapter 15, sections 
5 and 6. 

The Markov chain {£/, Ha"} is evidently irreducible (since p k > 
for k < m) and aperiodic. Therefore lim P°(k,K,N) exists. Since it is 

A'-. M 
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irreducible, the chain has either all transient, all recurrent null, or all 
recurrent non-null states. 

If a state (k,K) is transient or recurrent null, thenlim P(k,K,N) = 0. 

ff-*co 

Therefore, to show that all states are recurrent non-null it will suffice 
to show that for some state (k,K), lim P°(k,K,N) > 0. It will then 

Ar-oo 

follow that this is so for all states, and that Yl P°(MO = !• We look 

at the state (0,0): 

To see that lim P°(0,0,2V) > 0, we compare our system (with arbi- 

N-oo 

trary balking probabilities) to the special system for which m = 0, 
M — oo, H K = (always assuming the same input process). For this 
special system, write P{ Sjv° — K\ = P"(K,N), and take as initial 
condition: S(0+) — i -\- I . 

It is clear that for any system with M = oo, and with the same 
initial condition, 

P°(0,0,N) ^ P°(0,N), 

for each N, whence 

lim P°(0,(),N) ^ lim P°(0,N). 

JV-»oo N-*oo 

But it is known 3 that lim P°(0,N) > 0; thus 

N—oo 

lim P°(0,0,JV) = P°(0,0) > 

and all states are recurrent non-null. Hence, since the chain is also irre- 
ducible and aperiodic, it is ergodic. 

We now know also that a unique stationary distribution exists and 
that it coincides with the limiting distribution. From (6), we must 
have 

U°(k,R) = ZyAR)[U°(j,R) + V°(j,R- 1)]. (26) 

Denote by U 0,R the row-vector with components U°(k,R),0 ^ k ^ m. 
Then (20) may be written 



Thus, from (18), 



U°' R = V e ' R -V'*. (27) 
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Writing out (27) in components, and using (12), we obtain (24). 
We now prove (25). Denote by C (B) the event that the nth arrival 
overflows. Thus, 

b = rimP{C (n) j. 

n-»oo 

Now, 
P°(k,K) = Urn /Uv° = *, »/ - K) = lim P{£„ = k, S. = K | C (n) j 

.. P{g» - fc, E, = 7ij P |C (, ° I t, - fc, H„ = gj 
~ »J£ P{C<»>} 

But 

p{c {n) | &, = fc, a, = A! = p{c (n) i e, - fc} - q k . 

Therefore 

P»(/ c> ^) = »*(*.*> (28) 



and 



<z* E P(fc, JO 

c/°(fc,o) - E P°(fc,«) = — ^ 



K=0 



= g.p(fc) QJ}D< 

To complete the proof of Theorem 2, it remains to show that the 
binomial moments U°(k,R) uniquely determine the probabilities 
P°(k,K). This proof will be easier after we have discussed the stationary 
distribution at arrival moments, P(k,K), and we therefore defer it 

until then. 

It is sometimes convenient to work with the double binomial moments 

B°(r, R) = E (*) U°(k, R) 

C°(r, R) = E (*) V°(k,R). 
In terms of these, (24) and (25) of Theorem 2 become 
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(il 



B°(r,I?) = £ [f Jr (R) - gj T (R)}C°(j,R - 1) 

J=0 

(R = 1,2, •■■ .il/) 
B°(r,0) =rZ( UkPik). 

O A-=r V / 



2o: 



CiO.) 



Here we have used the following definitions: fi r (s) and gi r (s) are the 
/tli differences of <?>or(s) and ^or(s) : 



V*) = E(-D' -> ■ )**(«) 

>-o V/ 



yir(«) = E (-D'"M. -)*yr(a) 



where *,>(s) and¥; r (s) are defined, following Takacs [(T59), (TOO)], by 



(31) 
(32) 



*>(•) = £ f J PlcHjk(s) 



and must satisfy [(TGI) and (T62)] 

*»(s) 



and 



<M«) = 



**•(«) 



1 - V(S) 



1 



') +*y ir _x(s) 



(33) 

(34) 

(35) 
(36) 
(37) 



c r (s) c r _i(s) _v 

as well as the relations in /• implied by their definitions [see (T25)], 

*>(•) = g Q (A^ Pr )* y< (s). 

Examples of the application of the methods of this section will be 
found in Section V. 

III. THE STATE OF THE SYSTEM AT ARRIVAL INSTANTS 



3.1 Transient Behaviour 

The sequence {£„ , E„| is clearly a homogeneous Markov chain. We 
assume initial conditions £(0+) = ?', E(0+) = /, and seek the dis- 
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tribution P(k,K,n). We no longer restrict our attention to states (k,K) 
for which q k > 0, but consider all states (k,K), ^ k ^ m ^ «>, 
g if ^ M ^ oo. 

We shall prove the following: 

Theorem 8: The distribution P(k,K,n) is uniquely determined by the 
double binomial moments 

B<r*0-g£(J)©IWM i 

the latter are determined by 

B(rAl) -!■*.(*) (L) (38) 

(r = 0, 1, ...,i»;fi -0,1, ••■ ,M) 

B(r,R,n + 1) = ^ r+ «[B(r,B,n) + D(r - l,/2,n) 

+ C(r,R - l,n) - #(r,fl - l,n)] (39) 

(r = 0, 1, •• • , m; R = 0, 1, • • • , M; n = 1, 2, • • • ). 
Here 

C(r,R,n) = E E rifpJOxP^An) 

A-=r K=K W V>V 

Is 0(1 

tfWW = £ £ (J) (£) p/W(M» 

and all these quantities are denned to be zero if r < or ft < 0. 

Proof: If the arrival at t„ finds the system in the state (j,J), it may 
either get on the first-choice group, with probability p, , or balk on the 
first-choice group with probability q } ; in the latter case, it may get on 
the overflow group, with probability Gj , or balk there too, with prob- 
ability Hj . Thus the transition probabilities are given by 

p(j,J;k,K) = P{$«+, = K S»+i - K | £„ - j, E„ = /} 



Z>(r,K,n) = t t (*) (f ) p k P(KK,n) 



-fM*('tV (i -^© ir 



-"(l-a"")" 



-*w+i-* (40) 



OVERFLOW TRAFFIC 63 

Now 

P(k,K,n + 1) = Z £ p(j,J;k,K) P(j,J,n). (41) 

Substituting (40) in (41) and taking binomial moments with respect 
to both the first-choice and overflow groups, we obtain : 

+ *©[ <7 '( J « 1 ) + ff '©]} POy ' B) - 

Note that the quantity in braces in (42) is 






+*C-i)© + H0G-i)}- (43) 

Substituting (43) in (42), we obtain (39). 
For n = 1, we have 

PflkAl) = [clF(x) Q c -^(l - 0**(£) •"**(! - •-)"; 

taking binomial moments with respect to both trunk groups, we obtain 
(38). 

From the double, binomial moments, one obtains the probabilities 
P(k,K,n) by using: 

»r-A (r 



U(k,R,n) = g (-ir* f y fl(r,R,n) (44) 

and 

P(*A») = E (-I)* - * (|) I7(jyi,n). (45) 

Clearly P{k,K,n) = for fc + K ^ i : + I + n; it follows that the 
sums in (44) and (45) contain a finite number of terms for finite n, 
even if M = oo , and there are no problems about convergence. 

Equations (38) and (39) may be solved, in some cases, by means of 
generating functions; we give an example. 

Exam-pie 2: 

We consider the simplest possible case, in which 

q k = (k = 0, 1, ••• , ro - 1) 

q m = 1 
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In this case, 



M = «> 
H K = (K = 0, 1,2, ••• 

C(r,R,n) = B(r,R,n), 
E(r,R,n) = D(r,R,n), 



and 



D(r,R,n) = B(r,R,n) - OnB(m,R y n). 

Substituting (46), (47), and (48) in (39), we get 
B(r,R,n + 1) = <p r ^[B{r,R,n) + B(r - l,R,n) 

-( ™ j) B(m,R,n) + M B(m,fi - l,n)]. 



(46) 
(47) 

(48) 
(49) 



Let 



B(r,R,w) = E B(r,R,n)w". 



From (38) and (49): 



B(r,R,w) = 



W<Pr+R 



i\/I 



+ B(r - l,R,w) 



1 — W(p r+ R \_\rj \R 

~ (r - l) B{ - m > R > w) + (?) 5(m ' R " ^j ' 
The solution of (50) is 

, S\71s^.-f /Ay (i\ 
\f(m\ 1 ftWftW 



(50) 



B(r,R,w) = r r +j»(w) 



1 

rj +s -i(w] 



jSiVj/r*H«(«).T^/AT 



ft \j7 r y+ *-i(tt>) 



S S 



1 a=o W ft \j/ Tj+a-i (to) 



1 
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where we have defined 

IV(.)..II "EL., (r- 0,1,2, ...) 

iu(») ■ i. 

3.2 The Limiting Distribution P(k,K) 

Theorem 4-' The quantities P(k,K) = lim P(k,K,n) exist, are strictly 

n-»w 

positive, form a probability distribution independent of the initial state, 
and are uniquely determined by the double binomial moments B(r,R) = 

£ (j U(k,R), where U(k,R) = £ P£\p(k,K); the B(r,R) are given 

by 

[_j = r Lj+R }-r+\ Lj+R-l J (51) 

(r = 0, 1, ■■■,m;R = 0,1, -.M). 
Here 

m / 7 \ / T r \ 

C°(r,R) = T l (!:)(t)G K P (k,K). 
k-r \r/ \tc/ 

Proof: That the limits P(k,K) exist and are independent of the initial 
state again follows from the faet that the Markov chain {£„ , a,,} (n = 
1, 2, ■ • •) is irreducible and aperiodic. To show that the P(k,K) are 
strictly positive and form a probability distribution, we must show that 
there exists some state (k,K) such that P(k,K) > 0. This can be done 
by a method similar to that used in the proof of Theorem 2; we omit the 
argument. It follows that a unique stationary distribution exists and 
that it coincides with the limiting distribution. We express this stationary 
distribution in terms of the stationary distribution P°(k,K) in the 
following way : 

Consider the arrival which occurs at t„ (under equilibrium conditions) . 

It either overflows, with probability b, or does not, with probability 
(1-6). 

If it overflows, the probability that it encountered the state (j,J) is 

P°U,J). 

If it does not overflow, let us denote the probability that it encountered 
the state (j,J) by P*(j,J). 
We note that 

P(j,J) = bP°(j,J) + (1 - b)P*(j,J). (52) 
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Suppose that 0„+i = x. 

If the arrival at t„ encountered the state (j,J) and overflowed, the 
probability that the arrival at t„+i encounters the state (k,K) is: 



+ #Y^V*(1 - O*"] = «(*), say. 



If the arrival at t„ encountered the state (j,J) and did not overflow, 
the probability that the arrival at t„+i encounters the state (k,K) is: 

( j ^ l \~*0 - c-*y +1 - k (jX-* K (l - e-*V K = (!(x), say. (54) 

Taking account of both these possibilities, and removing the condition 
on 0„+i , 

P(k,K) -£f f dF{x)[bP°{j,J)a(x) + (1 - b)F*(j,JMx)l 
Using (52), 

P(KK) =ZE rfF(.X'){6P°(i,J)[ar(.T) -j9(*)l + P(j/)^(*)}. 

y=o j=o •'o 

Taking binomial moments with respect to both trunk groups, and 
using (53) and (54), 

Birfi) - *,» 1 5 (wov) [(j) (ft(' 7 + : ) + //, g)) 

+ 6[C°(r,fl - 1) - B°(r - l,fi)]}. 
The solution of (55) is 
B(r,B) _B(m,R) f & [g gOg) _ £ W-l) 1 (g6) 

Now note that, from (28), 

bB°(m,R) = 5(7/1,72). (57) 

Substituting (57) in (56), we obtain (51). 
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To complete the proof of Theorem 4, it remains to show that the 
double binomial moments B(r,R) uniquely determine the probabilities 
P(k,K). It is clear that the B(r,R) uniquely determine the U(k,R) 
through the equation 

U(k,R) = g (-iy~ k m B(r,R) (58) 

because m is finite. Thus we must show that 

P(MO = Z (" 1 )*~* (£) ^ <*,«) (59) 

when il/ is infinite; it will suffice to show that the series on the right 
converges absolutely. 
From (39) we have 

B(0,R) = T -^— [C(0,R - 1) - #(0, R - 1)]. (60) 

I — (PR 

Now, 

C(0,R) - E(0,R) = g E k (*) (^) q«G K P{k,K) ^ B(0,R). (61) 
Therefore, 



Now 



whence 



Thus 



#(0,ft) S -^—B(0,R - 1). (62) 

1 — <p/t 



lim?„ = limp(«) = F(0+) = 

R-»» 8-.O0 



<Pr 



lira — ^— = 0. 

K-oo 1 — ^>ft 



r B(Q t R) n /Ain 

b(o^i) =0 - (63) 



Equation (63) is sufficient to insure that 
converges. 
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Consider for simplicity the case m — 1. Then we have 

B(0,R) = U(0,R) + U(l,R). (64) 

At least one of the statements 

is p(M-D = ° (65) 

.. t/(l,fl) n / rr x 

SS c;(i,b - i) - ° (66) 

must be true, for if both failed to be true, then for some e > there 
would be terms for which 

17(0,22) e 



U(0,R - 1) 
U(1,R) 



> e 



U(l,R - 1) 
for arbitrarily large R; it would follow that for arbitrarily large R 

B(0,R) U(0,R) + U(l,R) > 

B(0,R-l) tf(0,fl - 1) + U(l,R - 1) 

which contradicts (63). 

Say (65) is true. Then the series 

converges; thus 

£©™»-£©«^ -£©"«*) 

converges, and this proves (59) for m = 1. The generalization to 
arbitrary m is straightforward. 

Corollary: We can now easily complete the proof of Theorem 2 by 
remarking that [using (28)] 

bU°(k,R) = &E 8 (fl) P W) 

-Sg)*P(M)sS(J)p(M)-l/(W) 

so that the series 
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p°(k,K) = £; (-i) M (|) pw*> 



converges absolutely, Q.E.D. 

We again defer examples to Section V. 



IV. THE STATE OF THE SYSTEM AT ANY TIME 

4.1 Transient Behaviour 
Let 

"' M //\ /k\ 

with Laplace transform 

0(r,R,s) = [ e""B(r,R,0dt. 

Jo 

Let M ik 1K {t) be the expected number of arrivals in (0,t] to encounter k 
trunks busy on the first-choice group and K on the overflow group, on 
the condition that £(0+) = i, E(0+) = /, with Laplace-Stieltjes 
transform 



*•"(«) = f e- SI dM ik ' K {x). 
Jo 



We also define several kinds of double binomial moments: 

m M /l.\ / J.'\ 

*tr 7 *(s) = Z E ( ' 1 ( H ] M.fc /A (s) 
k=r K=R V / V"-/ 

fc-r k=R V / W 

5V R ( S ) «tf (T) (T) P*<W*(s). 



A=r X=R V 



A' 



Theorem 5: 



**'*(*) = 



^(s + r + fi) 



1 - «,(« + r + R) 

[(;) (£) + *,>-!'*(*) + x,/- R - i ( S ) - f 1 -/- h - i c S ) 



., (67) 
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Proof: Consider a certain set of r first-choice trunks and a certain set 
of R overflow trunks. We shall call the union of these two sets an (r,R)- 
tuple of trunks, and if the r first-choice trunks and the R overflow trunks 
are all busy at time *, we shall say that this particular (r,R) -tuple of 
trunks is busy at time t. Thus, when the system is in the state (k,K), 

the number of busy (r,ft)-tuples is I M R )- Let us make the conven- 
tion that there is always one busy (0,0)-tuple. The expected number of 
busy (r./j?) -tuples at time / is evidently B(r,R,t). 

Let us now calculate the expected total number of encounters between 
arriving calls and busy (r,R) -tuples in the interval (0,/]. Denote this 
expectation by E ir IR (t). 

If the nth arrival occurs in (0,t], and if (£„ = k, a n = K), then the 

nth arrival encounters ( j ( p) bus y (r,/2) -tuples. Thus 

*•>'"« " S % £ (r) © C dP(T " = * f » = "> H " " Kh 

But 

Z P{r n ?Zu,Z n = k, S„ = K) = M ik IK {u). (68) 

n=l 

Therefore 

*-»-££©©»•"«> 

with Laplace-Stieltjes transform 

Cl / R ( S ) = *,/*(«). (69) 

But e iT IR (s) can be found in another way. If (£„ = k, S„ = K), then 
at time T n +, the system is in the state (k + l,K) with probability p k , 
the state (fc,X + 1) with probability q k G K , or the state (k,K) with 
probability q k H K . Thus the expected number of busy (r,R) -tuples at 
time t„+, under the stated condition, is 

and the expected number of busy (r,fl) -tuples created by the nth arrival, 
under the stated condition, is 
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Now the probability that the life of a busy (r,/2) -tuple will be longer 
than x is exp ( — (r + 7?).r). Thus the expected number of encounters 
between arriving calls and created (r,/?) -tuples in the interval (0,/j is: 

w m M r-t 

EE E / dp \ Tn ^ u. £, = fc,s„ = #} 

n=l *=r— 1 K=H— 1 ""0 

Jo 

where ^/(.r) is the expected number of arrivals in an interval of length 
x, when there was an arrival at the start of the interval. M(x) has 
Laplace-Stieltjes transform 



1 - r(s) ' 

Equation (70) is a convolution. Recalling (68), we see that (70) has 
Laplace-Stieltjes transform, 



JUK-0© 



+ (1 - p k )Gjc (fy ( R f X\ nJ K (s)n(s + r + R). 



(71) 



We must not forget the (r,R) -tuples which were busy initially; the 
expected number of encounters between arriving calls and these is 



— (H-Jt)H 



£) t (jpw. s «).-— - (;) g) { MMe 

n 
U(a + r + £). (72) 



with Laplace-Stieltjes transform 

(:)G 

Adding (71) and (72) we get 



IR f \ 

Ur is) = 



<p(s + r + R) 



1 - <p(8 + r + B) 

Now comparing (09) and (73), we obtain (07). 



(73) 
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Theorem 6: The distribution P(k,K,t) (t > 0) is determined by 
Proof: We have 

P(fc,K,0 = Q .-*(! - c-'y- (£) e-"(i - <r')'-*[i - P(01 

oo m M ft 

+ ZEE dP{f n = j, Sn - J, r. ^ w} 

n = l J— o y=o •'o 

•^ftO e "" u,i(i - e " ( '" u,)7+1 " fc © 



<e -(,-u)K (1 _ e -it-u)y-K + (J \ -lt-u* (75) 

\ / 

. (1 _ e -a-u) ) , + i-K + /// ^ c -(-)x 

. (1 _ e -^)^l|[l -i?0-it)]. 

This may be seen as follows: either no calls arrive in the interval {0,t], 
or the last call to arrive in that interval is the nth (n = 1,2, •■•), i.e. 
the nth call arrives at time u and no calls arrive in the interval (u,t\. If 
this call encounters the state (j,J) it may get on the first-choice group 
(probability pj), the overflow group (probability qfij), or neither 
(probability qjHj). Then enough calls must end in the interval (u,t] so 
that the state at time t is (k,K). 

From (75), and keeping in mind (68), 

B(rM = (;) © e-'^n - «oi + 1 1 ( «,,"« 

•Gii)}' 1 -^-^ 1 - 

and taking the Laplace transform, 
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(3(r,B,s) = 



1 - <p(s + r + g) 
s + r + 72 



A// 



r/\72 



+ **"(«) + *,>-l'*(s) 



+ X ir IlB - l (8) - F, 



r (S) . 



(76) 



From (76) and (67) we obtain (74). 

It remains to show that the double binomial moments uniquely de- 
termine the probabilities P(k,K,t). As in the proof of Theorem 4, it will 
suffice to show that for alll > 



r 73(0,72,0 

l ™ B(0,R-l,t) =°- 



From (67), for R > /, 



r* ' R r^ <- ^( s + ^) * i.R-lr \ 
®i0 W S — — * i0 ( s ). 

1 - <p(s + 72) 



(77) 



(78) 



But, for all s > 0, 



Therefore 



li m V (s + I{ ) « o 
«-» 1 — <p(s + 72) 



lllll 



= 0. 



Now from (74), 

/?(0,72,s) 1 - v (s + 72) *>(s + 72-1) 



(79) 



0(0,72 - l,s) p(s +7?) 1 - ^(s + 72 - 1) 

s + 72 - 1 ** /a (s) 



and so 



smcc 



lim _.?<***) _i im *."W =0 , 



«-«. /3(0,72 - l,s ) b-oo *« r -^ l («) 

v(s) 



s + 72 W+*{8) 

(80) 



lim 



= 1. 



g-oo <p(S — 1) 

From (80), and the inversion formula for the Laplace transform, the 
result (77) follows. 
Example 3: 



74 
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(ft = 0, 1, ••• ,ro- 1 



Consider the case 

q k - 

g. = 1 

H K = H, G K = G(G + H = 1) (K = 0, 1, 2, ■••)• 
This example may be of some practical interest. It represents a situa- 
tion in which some equipment, other than a free trunk, is needed to set 
up a connection on the overflow group. If this equipment is serving a 
large number of trunk groups, the chance of its being idle may be sub- 
stantially independent of the situation on the particular overflow group 
being considered here, and may be represented by a constant, G. 
In this case we have 

X ir IR (s) = Gft* n (a) 
Y ir in (s) = (»*."(•) 
and 



*,/"(*) = *,y' W (s) - 



(;)*.»(•>■ 



Equation (07) becomes 



*</*(«) - 



<p(s + r + B) 



3(3 + ^" w 



1 - *>(s 4- r + «) 

-( r r i )**» /a (a)+G(7)*^ s - 1 ( 8 )^ 

The solution of (81) is: 



(81) 



•*"(») - Cr+«(«) ^ 



fyM _J_ 

&\j)c» K (8) r/Ar^' 



e ; e ; 



A/m\ j~ _ ^ W '"° ^' GSC > +s - l(s) 



a/™\ i « w « w 



yW ! 

p&\j/G*C i+li -i(s) 



fcg W ^ Vj7 Cfd+a-iia) 



(82) 



OVERFLOW TRAFFIC 75 

The expression for p(r,R,s) can now be obtained from (82), using (74). 

4.2 The Limiting Distribution P*(k,K) 

Theorem 7: The quantities P*(k,K) exist, are strictly positive, form a 
■probability distribution, are independent of the initial state, and are uniquely 
determined by the double binomial moments 

the latter satisfy 

B*(r,lt) = -£— * ~ * r+R B{r,R), for r + R > (83) 

B*(0,0) = 1. 

Proof: To prove the existence, we consider the limit of (75) as t—> ». 
Clearly the first term goes to zero, and we have 

.1/ « »! 

P*(k,K) - lim s /E^V(«) 



e -«->*(l - e -<'-»>)'-*ll[i - F (t- 



(84) 



«)]. 



It follows from Smith's "fundamental theorem," 14 the assumption 
that F(x) is not a lattice distribution, and the fact that P(j,J) > for 
all j and J, that the limit in (84) exists and is given by 
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P*(k,K) = tt^ 1 ^- f du[l - F(u)] 
j=o j=o a Jq 

■L (' + l ) .-(l - O* 1 " 

•(£)«-*(l - "-")'-" + iJ^jo-^l - a-)'-" (85) 

It is clear from (85) that P*(k,K) > for all (MO, since the integrands 
are all strictly positive. (Note also that we have assumed a > 0.) The 
dependence on (i,I) has disappeared, and it is easy to show from (85) 
that 

ES P*(k,K) = 1. 

/.=0 K=0 

Thus B*(0,0) = 1. To show (83), we take a different tack: 

Consider any state (MO- Transitions into the state (MO are of four 
types: 

(fc - 1,JST) — > (MO (type a) 

(MC " 1)->(M0 (type 6) 

(fc + 1,K) -+ (MO (type c) 

(k,K + 1)->(M0 (typed). 

Transitions out of the state (MO arc also of four types: 

(k,K)-+(k - l,K) (type a') 

(k,K) -+(k,K - 1) (type/;') 

(MO -♦(*+ W (typec') 

(MO -> (fc,7C + 1) (typed'). 

Denote by N v (t) the expected number of transitions of type y in the 
interval (0,/]. 

If Ave consider the process only at times when the state (MO exists, 
transitions of type (a') form a Poisson process of density k, and transi- 
tions of type (b') form a Poisson process of density K. Thus, 

AUO = k ( P(k,K,t)dt (86a') 

■'n 
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N b >(t) = K I P(k,K,l)dl. (80//) 

•"0 



Similarly, 



N c (t) = (k + 1) f P(k + l,K,t)dt (86c) 

N d (0 = (K + l) [ P(k, K + 1, t)dt. (8Gd) 

Now {£„ = k, E n = if) is a recurrent event, with mean recurrence 
time [a/P(k,K)\ > 0. Thus, from the "elementary renewal theorem," 15 

Iim M «*«) = gfeg) 

<-.» t a 

But clearly, 

tf*(0 = qkG K M ik IK (t), 

so that 

lim N *'M = 9*G*P(W) = hG K P\k,K) 



Similarly, 

lim ^2 = GK-ibP°{k,K - 1) 
t-»w / a 

lim ^''W = P^MO = P(W) ~ l>P°(k,K) 
«-»« /a a 

lim 2s^ . P(A: ~ Lgj ~ bI "(k ~ hK) 

t—ao I a 



(86a") 

(866) 

(86c') 
(86a) 



We now notice that in any interval (0,/], the number of transitions out 
of the state (k,K) can differ from the number of transitions into the 
state (k,K) by at most 1. From this remark, and all the equations (86), 
it follows that 

(A- + K)P*{k,K) + aP(k,K) - abH K P (k,K) 

= ah[G R -iP {k,K - 1) - P°(k - 1,/v)] + aP(k - 1,K) (87) 

+ (fc + l)P*(k + l,K) + (K + l)P*{k,K + 1). 

Taking the double binomial moment of (87), one obtains 
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(r + R)B*(r,R) = a [b{v - 1,7?) - ( m + *) B(m,R) 

1) (88) 
+ 6[c'(r,fl - 1) - B°(r - l,R) + ( W + l ) *•<«■»*) JJ. 

We now note that, according to (51), 

a[B(r- 1,7?) - (™ + *) B(m,ft)] 

= aba+a _rt^-i;2^ii)] (so) 

- ab (™ + *) JT(m,fl). 

Putting (89) into (88), we obtain (83). 

It is now easy to see that the B*(r,R) determine the P*(k,K). For 
from (83) 

.. B*(0,R) ,. r+R-2 <p*-i B(0,R) 

22 B*(0,R- 1) " 22 r+fl vr B(0, R - 1) 

= 1 R^ B(0,fi-lj -°- 

Corollary: For Poisson input, P*(k,K) ** P(k,K). 

Proof: For Poisson input, F(x) = 1 - e "', < a < co ; a = 1/a. 

Thus 

/ \ a 



B *( rjR) = -J^L+l B (r,R) = B(r,R), 
r -\- K a. 

and since the double binomial moments determine the probabilities 
uniquely, the result follows. 

Examples will be found in the next section. 

V. EXAMPLES FOR THE STATIONARY PROCESS 

5.1 Categories of Examples 

In this section we will try to find the stationary binomial moments 
B°(r,R), B(r,R), and B*(r,R) for certain special cases, or categories 
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of cases. In the easiest cases we will succeed in finding explicit expressions 
for all these moments; in a harder case we will find explicit expressions 
only when R = lorR = 2;m the most complicated example (the ran- 
dom slip with overflow group, mentioned in Section I), the treatment is 
numerical, and only the results for the over-all blocking, B, are reported. 
If the first-choice group is full-access, the situation is particularly 
simple, since overflow can only occur if £„ = m; the vector equations 

(24) for U°(k,R) then become scalar, and B°(r,R) = ( m j U°(m,R). 

If the balking on the first-choice group is arbitrary, but the overflow 
group is infinite with no balking, or with constant balking probability, 
as in Example 3 above, some simplification occurs. For then, 

V°(k,R) = GU°(k,R) 

and hence (24) becomes a recurrence relation, although the quantities 
it relates are vectors. In such a case it is straightforward to find the first 
few moments of the distribution on the overflow group. 

In cases in which neither of the above simplifications occur, the form 
of the balking probabilities may still be such as to facilitate calculation; 
an example of this is the random slip with overflow group. 

5.2 Full-Access First-Choice Group 
We suppose 

2 , = (k -0,1, ••• ,m - 1) 

Qm = 1. 

E( (nations (24) reduce to the single equation 

U°(m,R) = p mm (R)V°(m, R - 1) (90) 

and from (13), 

y(R) is given by (1); it easily follows that 

± ( m ) -1— 

v ( m ) _L_ 
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Noting that, from the definitions, 



cm = Sa, 

IsR—l 



(91) becomes 



Hm 



.(A) = 



£H_L 

J=0 \J/ Cj+R-l 

h\j/c J+R 



(92) 



(93) 



We also know [from (25)] that 



0.(«JO) -£&£-!. 



(94) 



Example 4-' 

We now consider a slight generalization of the system considered by 
Brockmeyer (see Section I). Namely, let 

qt = (fc = 0, 1, ••• ,m - 1) 

?m = 1 

H K = H (K -0, 1, ••• ,M - 1) 

ff* - I- 
In this case we have 

nm,fi) - <?|V(m,fl) - (Y) ET(ro,M)]. 
Thus, from (90), 
U°(m,R) = » mm (R)G[u°(m,R - 1) - Q^ J C/°(m,M)] 

(22 = 1,2, ... ,M). 
The solution of (95) is 



U°(m,R) = 
Now, from (93), 



■ . 1 t( M r)\G J U^ m (Q)'] 

G R n»uQ)\ J ; f R ;Lr r 4=r- (*» 



SCOM"-"'! 



OVERFLOW TRAFFIC 



81 



n aw(q) - -, 



UDl 



Ml, 



J=0 \j / Cj +K 



Thus, 



Z ( m> ) - 

B°(r,R) = (™)Q* 



(R - 1,2, •••)• 



£i G' h\j)C, 



S(7)c£ 



j+j 



.1/ 



z^/z /A ' 



£o 6 v Uyj/c 



'i+J 



We notice [see (T54)] that 
Thus, from (51), 



P(m)' 



B(r,R) = G R C r+ n 



(?) 



E^/zr 



=• G J U\j/Cj 



m\ 1 



r fW_L £ fa) . I 

i-°r \J / Cj+R _ j=r+i \j / C'j+R-i 

£ (?) a 1 - , 

J=0 \J / (-J+R-1 . 



t r ) — 



(R = 1,2, ...,Jf). 



(97) 



(98) 



(99) 



5*(r,i2) follows from (83). 

When G = 1, (99) is the generalization to recurrent input of Brock- 
meyer's result, (4). It can indeed he verified that (99), for Poisson input 
and for 6=1, agrees with (4). 

For infinite full-access overflow group (.1/ = °o , f/ = 1 ), (99) becomes 



B(r,R) = C r+R < 



z r 



7/1 



i = r \J ' W+« J=r+1 \J / Cj+R— 



III 



1 



T/l 

E 



■///. 



1 



£n — i , 



(100) 



Equation (100) is the generalization to recurrent input of Kosten's 
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result, (2). Again it can be verified that (100), for Poisson input, agrees 
with (2). 

5.3 Constant- Balking Overflow Group 
We suppose that M = <*> 

G K = G (K = 0,1,2, •••)• 
Then (24) becomes 
U°(k,R) = q k G Z^(R)U°U,R ~ 1) (A = 1,2,3, •••)• (101) 

Example 5 

Suppose further that 

q k . = q (k = 0, 1, •■• ,m - 1) 

q m = 1. 

This might describe a system in which some auxiliary equipment is 
needed to set up a connection on the first-choice group, some other 
auxiliary equipment is needed to set up a connection on the overflow 
group, and the probability that the auxiliary equipment is idle is con- 
stant, but this probability is different for the two groups. This is a rather 
plausible system, except that the overflow group is infinite. 

We note that the blocking for such a system is 

B = E E q„H K P{k,K) = H[q + pP(m)]. 
It is easy to show by the methods of Ref. 2 that in this example 

B(rfi) - fC r %\L/j® (102) 

so that in particular 

P(m) = B(mfi) = m * - ■ 

Thus, 
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Instead of (101), wc use (29), which in our case becomes 



B°(r,R) =GZ [f jr (R) - gj r (R)]B°(j,R - 1) 

j=0 



(103) 



(R= 1,2, •••)• 



In this case we have, from (37), 

* /r (s) = pU Jr (s) - (7) **»(*)]• ( 104 ) 

We can solve (35), (36), and (104) to obtain 

Ms) = v m s 1 " (lS v) Pcha J LS (?) ^fcd 

" L§ v/^c^yJ'LS-. \i) p i ci-m])' 

It follows from (31) that 

'± ( m ) 1 
f r \k / p k C k (s) 1 



Ms) = p r C r (s) 



m 



Lti> \kj p*C k (s) 



1 p'Ci^(s) 



_ (105) 



1 



if I > r 



P l C^(s) 

if Z ^ r 

From (105), fi T (s) — gi r (.s) can easily be calculated by observing 
that in this example 

fir(s) - g lr (s) = qf lr (s) + p OfjfUa). 
Then, from (103) one obtains 



B°(r,R) = G 



[>(:)+ ««wo£(:U£ 



P k C k (R) 



£»\k) p"C k (R) 
,_o p'Cj-iilx) i=r+1 p'Cj^i{R) 



(10G) 
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Noting that, from (30) and (102), 



B°(r,0) = 



qB(rfi) + V [ ™ ) B(m,0) 



in 



q + pB(m,0) 



(107) 



m 



1 



+ P 



K kJ p k C k 

we can use (106) to find B°(r,l), B°(r,2), etc., and in particular, the 
first and second moments of the distribution on the overflow group only, 
at overflow instants, £°(0,1), B°(0,2). The formulas are long; we quote 
only : 



£°(0,1) = G< 



<A + 






(108) 



( 



5.4 Other Cases 

Once B°(r,R) is known, it is straightforward to determine B{r,R) 
and B*(r,R), using (51) and (83) respectively. [If B°(r,R) is known, 
C°(r,R) can be determined, for use in (51), from the relation, which 
follows from their definitions: 

r(r,R) - £ Q (A r " t O R )B (r,J) ] (109) 

see (T45).] The problem is thus to determine B°(r,R), from (29) and 
(30), or equivalents to determine U°(k,R) from (24) and (25). We 
consider the latter method. 

To use (24) and (25), one must first of all determine Mj*W for all 
relevant j, k, and R [say, from (T70)], as well as P(k) [say, from (T44) 
and (T45)]. Then the V°(k,R) must be expressed in terms of the 
U°{k,R); in general V°(k,R) can be expressed in terms of the U°(k,J), 
with J ^ R, by a relation analogous to (109): 

V(k,R) = Z (jQ (A J - R G R )U°(k,J). (110) 

When (110) is substituted in (24), one obtains a set of simultaneous 
equations for the U°(k,R). Equation (25) serves as a boundary condi- 
tion. If M is finite, (24) can be used to express U°(k,M - 1), 
U°(k,M — 2), ••• , U°(k,0) successively in terms of U°(k,M), and 
(25) can then be used to determine U°(k,M ). 
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When the U°(k,R) are known, one finds the B°(r,R) by taking bi- 
nomial moments, and then the B(r,R) from (51). The probabilities 
P(k,K) then follow by inverting the binomial moments, and the over-all 
blocking is determined by 

m M 

B = ZZ<lJlKP(k,K). 

i=0 K=0 

Example 6 

We consider the system described in Section I 



/// 



-0/ 


'(:) 


(k = 0, 1, 


II K = 




(K = 0, 1, 


H M = 1. 







M - 1 



The IBM 7090 computer at Murray Hill was programmed to find the 
blocking probability B for certain values of the parameters, namely: 

m + M = 10 

y + M = 6. 

The calculations were carried out for two kinds of input: 

(i) Poisson 

(ii) That sort of recurrent input which is itself the overflow from a 
group of m trunks to which a Poisson stream of calls (with negative- 
exponential holding times) of mean intensity a is submitted. Note that, 
since Poisson traffic is completely characterized by one parameter (its 
mean, in our case a ), this sort of recurrent input is completely charac- 
terized by two parameters (an and m ). 

Note also that this program allows one to calculate B for certain more 
complicated trunking arrangements, in the case of Poisson input, e.g., 
2 common trunks overflowing to a random slip of 3 on 7 which in turn 
overflows to 1 common trunk. (This arrangement also involves a total 
of 10 trunks and f> crosspoints per line. 1 

The results (blocking probability B as a function of input traffic a) 
are shown in Tables I and II and Fig. 1 . The cases treated were ?n = 
(Poisson input, a — a ) and mo = 2, in which case, of course, 

a = t/ V + a ° + T 

7 was given the values 2,3,4,5,0. ( Note that if 7 = 6, then M = 0; there 
is no overflow group.) 
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I 2 3 4 5 6 

SUBMITTED TRAFFIC, a, IN CALL- HOURS 

Fig. 1 — Blocking, D, vs submitted traffic, a. 

Before commenting on the results, we mention parenthetically several 
special features introduced into the calculation by the special form of the 
balking probabilities and by the kind of input process considered in this 
example. First, as to finding the P(k): (T44) and (T45) read, in our 
notation 



B(r,0) = 



1 - 



D(r - 1,0) 



D(rfi) = g (i\ (A*»BO0>. 



In the present example, 

<Pr m «o C r m "(flo) 

i - & ' r a + r°(oo) 



(r = 1,2, ••■) 



(HI) 
(112) 

(113) 
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and 

(,' ) 

A ^p r= _y7/ (j--r+l,f + 2 f --). (H4) 

(;) 

Also, since the overflow group is full-access (although finite), the 
relation (110) becomes 

V°(k,m = U°(k,R) - (¥\ U°(k,M). (115) 

In Tables I and II and Fig. 1 , we have used the notation y/m + M 
to describe a random-slip configuration in which each line has access to 
7 out of the m first-choice trunks and all the overflow trunks, except that 
the case 7 = G, m = 10, M = is referred to as 6/10. The curves in 
Fig. 1 have been drawn, to avoid crowding, only for 4/8 + 2 and 6/10. 

The following conclusions can be drawn from these results: 

(i) The blocking is higher, for the same mean traffic, when m = 2 
than when m = 0. This is consistent with the intuitive notion that 
overflow traffic is "peaky". 

(n) In a practical range of blocking (B = 0.001 or 0.01), 4/8 + 2 
is the "best" arrangement and 6/10 is the "worst" of those considered, 
from the point of view of the traffic capacity of the system for a fixed 
blocking probability. It can be seen from the curves that if one wanted 
an arrangement using 6 crosspoints per line and 10 trunks, one would 
gain about 8 per cent (for m = 2) or 6 per cent (for m = 0) in traffic 
capacity at B = 0.01, by using the arrangement 4/8 + 2 instead of 
6/10. At a blocking probability B = 0.001, these gains would be about 
16 and 11 per cent respectively. Such increases in traffic capacity are 
not negligible; they seem to be larger for peaky traffic than for Poisson 
traffic. 

(iii) For higher blockings ("overload" conditions), the advantage of 
4/8 + 2 relative to 6/10 diminishes. 

A study for a practical case would involve calculations of the block- 
ing for other values of 7 + M, a knowledge of the relative costs of trunks 
and crosspoints, and of course many other considerations, such as the 
relative costs of building and controlling 4/8 + 2 and 6/10 switches. 
Also, in such a study, one would want to keep in mind the approxima- 
tions implicit in the model used in this paper. For example: 

(i) In reality, blocked calls may wait or be resubmitted. 

(n) In reality, the number of traffic sources (lines) is finite, so that 
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the arrival process after any instant is dependent on the number of 

trunks busy at that instant; thus the input is not, in reality, recurrent. 

(in) As a further result of the finitcness of the number of lines, the 

complete set of ( 1 access patterns required for a perfect random slip 

probably could not be used, and even if it could, equal traffic would not 
be submitted to each access pattern (so that the blocking experienced 
by different subscribers would be different). 
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